HIGHER CONSERVATION LAWS FOR THE QUANTUM 
NON-LINEAR SCHRODINGER EQUATION 

B. Davies and V. E. KorepiiJ^ 

School of Mathematical Sciences, Australian National University, 
GPO Box 4, Canberra, ACT 2601, Australia 
(Dated: August 1989) 

Abstract 

We construct explicit forms for two non-trivial conservation laws of the quantum non-linear 
Schrodinger equation and show that they have the correct quasi-classical limit. For the second 
quantised form cannot be obtained by normal ordering of the classical conserved quantity H'^. We 
show that the Quantum Inverse Scattering Method also gives the correct higher Hamiltonians H^, 
and H/^. The surprising result is that the expansion of fundamental integrals of motion such as 
j4(A), in inverse powers of A, cannot be recovered by normal ordering of the classical expansion. 
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I. INTRODUCTION 



The quantum inverse scattering method (QISM) has its origins in attempts to extend 
the classical inverse scattering method for integrable non-linear systems [l-^ to interacting 
quantum fields j4-7|. Initial investigations, which were mainly concerned with semi-classical 
quantisation, were soon followed by a scheme for the exact quantisation of the non-linear 
Schrodinger equation jsj. In the last ten years there has been rapid progress on the QISM, 
with many published papers on the subject. Three reviews to which the reader may refer 



are refs. |9l-lll| . Most papers have been concerned with the development of this new branch 
of mathematical phys ics^ although there have been some whose concern is with the math- 
ematical foundations 
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13| . or which have raised objections particularly for the 



example of the quantum non-linear Schrodinger equation (QNLS). The most serious of the 
proposed difficulties relate to the higher conservation laws: it is the main concern of this 
paper to resolve these questions, at least for the QNLS. 

Let us define the model and recall some well known results 16[]. We consider the quantum 



non-linear Schredinger equation in 1 + 1 space-time dimensions. In second-quantised form, 
the Hamiltonian is given as 



2 = dx + c\l>'^(x)\l/^(x)\l/(x)\l>(x)}. 



(LI) 



Here \l/(a;,t) is a Bose field satisfying the canonical commutation relations 



[*(x), ^\y)] = S{x - y), [^{x), ^{y)] = 0. 



(1.2) 



Operators for the number of particles Q and total momentum P are given by 



(1.3) 



They are integrals of motion: [if2,-P] = [H2:Q] = [PiQ] = 0. For the repulsive case 
c > 0, the only case considered herein, a complete set of eigenfunctions of the operators are 
well-known: viz 



Ai,...,A^) = (iV!)-i/2 / d''xxN{xi,...,Xn\Xi,...,XN)'^\xi)---¥{xN)\0), (1.4) 



where the exphcit formula for the functions xn is 



Xn{,Xi, . . . ,Xn|Ai, . . . , \n) 

N 



^(-1)^ JJ{Apj - Apfc -zcsgn(xj - Xfc)}exp 
p j>k 



(1.5) 



n=l 

Here P runs over the permutations of (1, ... , A^) and sgn(x) is the sign of x. 

In addition to the second quantised form fll.ip there is a formulation of the QNLS as a 
non-relativistic many-body problem, in terms of partial differential operators and boundary 
conditions. In the A^-particle sector, the functions xn arise as eigenfunctions of the following 
differential operator H2: 

^ ^2 

^2 = -E^ + 2c E ^i^-^'^)- (1-6) 

j=l i N>j>k>l 

The delta function interaction may be replaced by boundary conditions at Xj = x^ and this 
will sometimes be done in the following. The functions xn are also eigenfunctions of the 
momentum differential operator P = —iHi. 

j = l J 



The eigenvalues of these two operators are given by 



N 



As we have mentioned, there are claims . Il5| that the higher conservation laws obtained 
from the QISM are in conflict with those which may be found directly from the above 
solutions, even for the next two operators and H4. We shall show that this is not the 
case. In the differential equation formulation we have only to correct an error in ref. fl^ . 
This we do in sections IlIII and UTTt where we construct the operators H3 and H4. They have 
the same eigenfunctions (II. SP with the eigenvalues ^ A| and ^ Aj, respectively. We have 
broken the calculation into two sections: in section [ITl we deal only with the N = 2 and 
N = 3 sectors where the calculations are quite elementary and already reveal the flaw in 
ref. [1^ . In section IIIII we give the forms in the general A^-particle sector. 

One must be most careful when writing down the corresponding conservation laws by 
means of quantum Bose fields because the individual terms in the formal expression (as a 



sum) involves irregular (undefined) operators, even though the total expression is very well 
defined. Expressed in the language of second quantisation, the results of sections [TTl and UTTl 
for may be summarised in the expression 

H, = jdx {^\x)<iI,Ux) - {3c/2)^\x)\^{xf),}. (1.9) 

There is no corresponding expression for H4. In particular. 



H,^ J dx {^L^..{X) + 2c(^t(^)2)^(^(^)2)^ 

+ c^^(x)2*^(x)2 + c^iixf'^ixf + 2c^\xf'^{x)^}. (I.IO) 

One may write down a formal second quantised form for H4 by replacing half of the three- 
particle interaction term 2c'$''{x)^'${x)^ by c\E'''"(x)^^(x)^^(x)\E'(a;)^, which is not normally 
ordered. The application of this symbol to a Fock space state such as (II. 4 p gives rise to 
the meaningless product 5^(xi — X2) of generalised functions, so it is hardly a useful mod- 
ification. At the same (formal) level of discussion, notice that the numerical coefficient of 



s us that any (formal) 
We shall see in sec- 



this interaction term is 2, the same as in the classical case. This tel 
quasi-classical limit will be correct, contrary to the claim of ref. jl4 
tions HI] and IIIII that H-^ and H4 are properly represented in terms of irreducible parts J3 
and J4 together with multinomials in the lower conserved operators. From H4 on there is 
no way of regrouping the formulae to give a second quantised form as the one dimensional 
integral of a density. Thus we cannot obtain Hn{n > 4) by normal ordering of the classical 



expressions which are integrals of such densities. In this respect we agree with Gutkin 
pages 112-114. 



The connection o 
is claimed in ref. 



the QISM to higher conservation laws is a more technical problem. It 
15| that the QISM fails to generate the correct conservation law even for 
H3. We shall show that this is not so: the difficulty lies in the asymptotic analysis. The 
anchor point of the QISM derivation of higher conservation laws is the fact that the trace of 
;he monodromy operator r(A), where A is the spectral parameter, gives a commuting family 

HI: 

[r(A),r(/i)]=0. (I.ll) 

This is true for both the lattice version and the continuous limit of the QNLS in a finite 
box: it has its analogue also for an infinite box. (Recall that r(A) is the transfer matrix in 



statistical mechanics.) In section HVl we discuss the quantum trace identities for the lattice 
QNLS. We show that the higher terms in the A~" asymptotic expansion are not given by 
normal ordering of the corresponding classical expressions. Our calculations show that the 



15|, are due to the neglect 



discrepancies in the asymptotic expansion, as reported in ref. 
of quantum corrections (contributions from operator reordering). We derive the large A 
expansion 



A— >— ioo 

where the commuting constants Aq to ^3 are 

^0 = —icHQ = ~icQ 
^2 



1 + \~'Ao + A-'Ai + X-^A2 + A-M3 + 0{y 



Ai = -cHi 



-HoiHo 



mi 



IC 



A2 = -icH2 + ic\Ho 
A3 = CH3 - "^Hl + c^d - H^)H2 
+ ^(iJo-l)(i^o-2)ifi + 



Ho{Ho-l){Ho-2) 
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Ho{Ho-l){Ho-2){Ho-3). 



(1.12) 

(1.13) 
(1.14) 
(1.15) 

(1.16) 



These results differ from those given in ref. 



15| in a number of respects. A minor difference 
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13| on 



is that our definition of A(X) follows the usual one of all the preceding literature 
the QNLS, whereas ref. 15|] interchanges the meaning of A{X) and its (Hermitean) adjoint 
A''" (A). More important are the differences in the operator ^3. First ^3 is not the normal 
ordered version of the correspond classical quantity. Second, we have corrected a numerical 



error in the coefficient of the third term of ^43: in [15| this is given as (2 — Hq). This 
correction is important because it makes it evident that the difference between the correct 
result and the normal ordering recipe stems from difficulties with asymptotics, rather than 
a fundamental flaw in the QISM. 

Finding the correct asymptotic expansion directly for the continuous model in an infinite 
box is an extremely tricky business indeed. The individual operators, A{X), v4"''(A), which 
appear on the diagonal of the monodromy matrix are themselves constants of the motion 



13| . One needs to expand the operator A{X) in inverse powers of A. The fact that the 



higher conservation laws are not normally ordered is equivalent to the surprising result 

5^A-'^-iA„^5^A-'-i:A^, (1.17) 



n>0 



n>0 



even though A{X) is correctly defined as :A{XY^:, the normal ordering of the classical quan- 
tity. The QISM has never depended, for its validity, on the normal ordering recipe. However 
it was difficult to see how this recipe could be broken for the continuous QNLS, even though 
there is no a priori reason to require it. In section |V] we will show that, when we con- 
sider A{X) as an integral operation in Fock space and take the asymptotic decomposition in 
inverse powers of A, the expansion is non-uniform and this leads to the breakdown of the 
formal expansion for n > 3. 



II. TWO AND THREE PARTICLE SECTORS 



In this section we first discuss two particle wave functions. They are given by 

X2(a;i,X2|Ai, A2) = {A2 - Ai - icsgn(x2 - xi))} exp{ixiAi + ix2X2} 

+ {A2 - Ai + icsgn{x2 - xi))}exp{2a;iA2 + ix2Xi}. (II. 1) 

This is a continuous symmetric function of Xi and X2', it is an eigenfunction of Hi 



Hi = d/dxi + d/dx2, 
H1X2 = «(Ai + A2)X2- 



(11.2) 
(11.3) 



It is also an eigenfunction of the Hamiltonian (II. ip 



Ho 



dx\ dx\ 



2c5{xi — X2) 



H2X2 = (A? + A^)x2. 



(11.4) 
(11.5) 



The operator (III. 411 is the free Hamiltonian except at the boundary xi = X2 where the 
interaction is equivalent to the following boundary condition 



CX2 + 



d 



d 



X2 



0. 



(H.6) 



dxi 8x2 

In the two particle sector these are the only two independent conservation laws: all higher 
conserved operators are generated by Hi and H2. It is easy to construct them. For conve- 
nience, introduce an operator J2 as 



J2 



dxidx2 



+ c5{Xi-X2) = UHI + H2). 



(11.7) 



The wave function X2 is an eigenfunction of J2 



-^2X2 = -(AiA2)x2- 

Now let us construct the operator ifs, with eigenvalues equal to i^(Af + A 



(11.8) 



(11.9) 



From (jll.sp . ( III.SP and ( ]11.7|) it follows that the Bethe wave function X2 is an eigenfunction 
of i?3, moreover an elementary calculation shows that it has the expected eigenvalue 



H3X2 = + A^)X2 



An explicit formula for is 



+ TT^ — OC0[Xi — X2 



d 



dxi 8x2 



(11.10) 



(11.11) 



and this shows that it coincides with the third conservation law constructed in refs. HQ. 



Now let us construct the operator H4 with eigenvalues equal to (A^ + Ag) 

= Hl + 2Jl-AHlJ2. 



(11.12) 



We emphasise that, by its very construction in terms of the lower conserved operators, 
is well defined as an operator. From ( jll.SP and (jll.SP it follows that 



H,X2 = (At + A^)X2. 



(11.13) 



So we have constructed the fourth conservation law for the QNLS in the N = 2 sector. It 

n 

does not coincide with the fourth conservation law constructed in [IJ]. Let us denote the 
latter by G4, it is given in [l^ as 



QA QA / g2 Q2 Q2 

Ga = TT^ + TTT - 2C(5(Xi - X2) I TT^ + TT^ + 



dxf dx\ 



dx\ dx\ dxidx2 

52 92 



2c 



+ 



+ 



92 



dxi dx\ dxidx2 



An elemenetary (formal) calculation shows that 

G4 = Hi-2c^S'^{xi-X2). 



S{xi-X2). (11.14) 



(11.15) 



This equation has only a formal significance, since 6'^{xi — X2) is an undefined product. So 
G4 is irregular (undefined) because is regular. This is clear even from (111.14^ because 



{d'^ /dx\)x2 ~ —X2)x2, so the product S{xi — X2){d'^ / dx\) is not defined as an operator. 
G4 does not commute with the Hamiltonian H2 and the Bethe eigenfunction X2 is not an 
eigenfunction of 6*4. We can already see the reason for (ll.iop . The operator (111.140 is 
precisely what we recover by using the normal ordered form of fll.lOp in the two-particle 
sector: for = 2 there can be no three particle interaction as contained in the term 
On the other hand, if we use the normal ordered symbols (II. ip and (II. 3P in 
eqs. (III. 71) and (111.12^ . we find that normal ordering cannot be carried out to rearrange the 



formula as the one-dimensional integral of a single Hamiltonian density. This is the meaning 
of the difference between G4 and H^^. 

We have constructed two non-trivial conservation laws in the two-particle sector, and 
used them to check the consistency of three operators H^, H4, G4. Now we discuss the three 
particle sector. The Bethe wave function xs is given by 



X3(a;i,X2,X3|Ai, A2, A3) = ^(-1)^ Y\ {\pj-\pk-icsgii{xj-Xk)} exp 

P 3>j>fe>l 



■ 3 

^y^XnApn 
n=l 

(11.16) 



It is a continuous symmetric function of Xi, X2, x^. It is an eigenfunction of the operators 
ifi, defined in the three-particle sector as 



3 o 3 02 



j=l ^ j=l J 3>j>fc>l 

The operator J2 = {Hf + -f^2)/2 now has the representation 

3>j>fc>l ^ J ^ 

It is well known that Xz is an eigenfunction of Hx and H21 and that for the Hamiltonian this 
property is equivalent to the boundary condition 



d 9 , 

CX3 + ^ ^ X3 

axj oxjj^i 



= 0. (11.19) 



It follows from our construction that it is a properly defined operator in the three particle 
sector and that xs is also an eigenfunction of J2: 



-^2X3 

^3>i>fc> 



E ^,^k]x3- (11.20) 
:7>fc>i / 



The third conserved Hamiltonian H-^ was constructed in Il5| correctly as the differential 
operator 

«^ = Egg>-^'= E ^f^^-^-') • (■■■21) 

j=l ^^j 3>j>k>l ^^^J ^^^^ 

Notice that the differentiation in the second term acts in a direction orthogonal to the 
argument of the delta function. This is important, since the Bethe wave function has 
discontinuous derivatives at the boundaries Xj = Xk- We rewrite in the form 

H3 = Hi - 3H^J2 + 3J3. (11.22) 

Here we have introduced the new operator J3 as 

d d d 

^3 = ^ — 5 — 5 ^ c——5{x2 - X3) + C-— 5(x3 - xi) + c-—5{xi - X2). (11.23) 

OX1OX2OXS 0x1 0x2 OXs 

From the construction of J3, we must have the equality 

^3X3 = ^'(AlA2A3)X3, (11.24) 
and this is equivalent to boundary conditions of the form 



_d_ 

dx^ 



d d 

CX3 + 1 ^ ^ ) ^3 

0X1 0x2 



= 0. (11.25) 

X2=xi+0 

which follow immediately from ( III. 191) . This shows directly that ( III. 241) is valid and that 



(III. 211) is the correct form for H3. In the three particle sector, all higher conservation laws 
are generated by Hi, H2 and -^3. We shall write them as functions of Hi, J2 and J3. Let us 
construct the fourth conservation law for an operator H4^: 

H^ = H^ + 2Ji - AHIJ2 + 4H1J3. (11.26) 

From our previous results it follows immediately that this a properly defined operator and 
that its action on the Bethe eigenstates is 

i^4X3 = (At + A^ + A^)X3. (11.27) 

As in the two particle sector, it does not coincide with the fourth operator G4 published in 
ref. [1^. Elementary (formal) manipulations show that in the three particle sector, 

Gi = Hi-2c^ ^ 5'^{xj-Xk) + Qc^5{xi-X2)5{x2-xz). (11.28) 

3>j>fc>l 

So 6*4 is not an integral of motion, it is not even defined for N = 2 and 3. Again the 
differences between G4 and H4 may be (formally) viewed as an ordering problem but it is 
not a profitable approach. 
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III. MANY PARTICLE SECTOR 



Recall the formula f ll.SP for the particle Bethe eigenstates XNixi\Xj). They are contin- 
uous symmetric functions oi Xi, . . . ,xn and Ai, . . . , A^v, and also eigenstates of the operators 
Hi, H2, defined in (II. 6 P and (II. 7p . This latter fact is equivalent to boundary conditions of 
the form (111.61) . The third operator H3 in the sequence of conserved quantities is given in 
Q, 3] as 

j=l J N>j>k>l \ J 

Let us write it in the form 

= Hi - 3H1J2 + 3J3, (III.2) 

as in section [ITl Now we have extended the definition of the operators J2 and J3 to the N 
particle sector as 

N>j>k>l \^-^3^-^k J 

f d d d \ 

N>j>k>l>l ^ J J ' 

(111.4) 

From ( II. SP it follows that 

■hXN = - 



^j^k]xN. (III.5) 

>i>fc>i / 



To prove that xn is also an eigenfunction of H^ we may first prove that it is an eigenfunction 
of J3: 

J3XN = i'^ I 

\N>j: 

and this is readily reduced to the following boundary condition, analogous to (111.250 . 



Y ^^^kXijXN, (III.6) 

j>k>l>l J 



y — 



d d 
cXn + I 7^ 7^ I Xn 



dxj dxj^ij Jxj+i=xj+o 



0. (III.7) 



So we have proved that 

HsXN = t'(j2x'i^^N, (III.8) 
and H3 is the correct operator for the third conservation law. 
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Let us construct the fourth conservation law in a similar manner. We commence with 
the definition 

Hi = Hi + 2Jl - AH^J2 + 4H1J3 - 4J4. (III.9) 
Here we have introduced the operator 



Ja 



N>j>k>l>m>l \ J ' 



+ c 



+ 



92 



+ 



dxjdxi 



6{Xk - Xm) + 



dxjdxk 
92 



5{Xi - Xm) 



6{xk - xi) + 



dxkdxi 



5{xj Xm) 



dxkdxm 



6{xj - xi) + 



92 



dxidxj 



-6{xj - Xk) 



6{Xj - Xk)S{xi - Xm) 



+ 6{Xj - Xl)6{Xk - Xm) + S{Xj - Xm)S{Xk - Xl) 



To prove that 



is equivalent to showing that 



N 



(III.IO) 



(III.ll) 



JaXn 



\N>a>b>c>d>l 



^a\^c^d XN, 



and this reduces to the following boundary conditions: 



TV 

N>j>k>3 



cXn + 



d 



d 



dxi dx2 



Xn 



0. 



(III.12) 



(III.13) 



Again these are valid, and follow immediately from the simpler boundary condition (III.19P 
because the differentiations are always in an orthogonal direction to the planes on which the 
delta functions have their support. In this way we have constructed the fourth conserved 
quantity H4,. It does not coincide with the operator of ref. [l^ (eq. 2.20), which we shall 
call G4. We repeat the formula here: 



2<= E 



ISC^ ^ S{Xj - Xk)S{Xk - Xl) 
J N>j>k>l>l 

Q2 Q2 Q2 



.dx"] dxi dxjdxk 

N>j>k>l^ 



Q2 Q2 Q2 

6{Xj - Xk) + S{Xj - Xk) { TT^ + TT^ + 



dxj dxf, dxjdxk. 



(III.14) 
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Comparing this with H4 we see that 



G4 = -^4 — 2c^ 6'^{xj — Xk) + 6c^ S{xj — Xk)S{xk — xi). (III. 15) 

N>j>k>l N>j>k>l>l 

So we have constructed ifs and H/^. We would hke to discuss their second quantised 



form and compare these with the corresponding classical conserved observables Hf and Hf. 
The latter are given in eqs. flIV.7l) . The operators J2, J3, J 4 are well defined as differential 
operators, and also in second quantised form. The square of J2 is also a well defined operator, 
but to carry out the desired comparison we need to write it as the sum of irregular terms 
(the whole sum is well defined). One of the terms involves (formally) the square of a delta 
function, another 5{xi — X2){d'^ /dxf), and the singularities cancel. So we can write H-^ and 
H4 in the form of a multi-dimensional integral of the field operators and their derivatives: 
this will be well defined but the individual terms in a rearranged one-dimensional expression 



may not be so well defined. With this caveat, for we recover the result of ref. IJ]: we 
could find the same result by applying the recipe of normal orderin g to Hf-see eq. 3.16. 



In the case of H4 however, we do not reproduce the result of ref. which we called 
G4 above. is not the normal ordering of H^^, although its (formal) quasi-classical limit 
is correct. This is an important correction to the result claimed in ref. |l^. As we have 
already observed, the QISM does not depend on the recipe of normal ordering, and in fact 
it is broken at the operator H^. However, it is required that the quasi-classical limit of an 
integrable quantum theory (when problems of ordering go away) should be correct. 



IV. QUANTUM TRACE IDENTITIES 

First let us discuss the trace identities as they were constructed in refs. We 
commence with the classical U operator 

U{x\X) = -^ + '-^a3 + Q{x), 
ax 2 



Qix) 



iy/cip^{x 

-i^/cip{x) 

The transition matrix T{x,y\X) is defined as the solution of an initial value problem: 



(IV. 1) 



U{x\X)T{x,y\X)=0, T{y,y\X) 



1 
1 



(IV2) 
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Suppose we impose periodic boundary conditions in a box of length L. Then the monodromy 
matrix T(A) and its trace r(A) are defined as 



T(A) = T(L,0|A), r(A)=trT(A). 



(IV.3) 



We are interested in the decomposition of exp(— zAL/2)r(A) in inverse powers of A, as 
A — —too. Let us make a gauge transformation which diagonahses T{x,y\X): 



T{x,y\\) = Vix\\)D{x,y\\)V~\y\\). 



(IV.4) 



Here D is a diagonal matrix while \^(x|A) and V ^{y\X) depend only on one space variable. 
l^(x|A) can be represented in the form 



V{x\X) 



(IV.5) 



1 /(x,A) 
/>,A) 1 

We note that /(x,A) — ?■ in the limit x — ?■ —ioo. If we substitute (IIV.4p and fllV.Sp into 
(I1V.2I) we get an equation for f{x,X) and an equation for D{x,y\X), from which it is very 
easy to get the 1/A decomposition of f{x, A) and for r(A). Here we simply quote the results 

.2 



'iXL/2 



r(A) 



+ 



A— >— ioo 
1 



A 



rcl\2 



A3 



^cHf + ^^H^'Ht + '-^{H^') 



cl\3 



1 

A4 



cU2 



c'H^'Hf + "^{H^'fHt + Y^iH^y) + ^(l/A^). (IV.6) 



We have introduced the following notation so as to keep the classical and quantum cases 
parallel 



^0^ 



m 



Hi 



dx ^lj{x)ip{x), 
dxip{x)ij^{x), 

dx {ipx{x)i'x{x) + cipixY'ipix)'^}, 
dx {4j{x)tlJxxx{x) - {3c/2)tP{x)\4j{x)% 

dx {i'xxix)'4)xx{x) + 2c{ilj{xY)x{'4'ixy)x 
+ cijixfipxixY + ctpxixYtpixy + 2ctlj{xft/j{xf}. 



(IV.7) 



-} 
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It is well known that the asymptotic expansion of the logarithm takes a simple form: in our 
present notation it is 



log(e-^^^/V(A))] ^ 1 - ^< - ^^Hf - '-^^Hf + ^^Hf + 0(1/A^). (IV.8) 
Now let us discuss the quantum case. One can repeat similar calculations in the quantum 



case to those above, for example this is done in ref . 2l| . Now one has the problem of non- 
commutation together with the fact that some of the necessary operations involve repeated 
differentiation of formally defined objects. The expansion up to 0(A^^) is 



J A— >— ioo A A \ Z 



+ 1 : [-tcHt + tc'H^'Ht - '-^{H^'Y^ : +0(1/A^). (IV.9) 

It is well known that everything is correct to this point. Performing the indicated normal 
ordering on the coefficients found thus far (with the help of the canonical commutation 
relations) we find that the first three commuting constants are given by eqs. ( II.13P - (lI.15p . 



Now let us take the next coefficient from the classical expansion flIV.6P and apply the same 
normal ordering prescription: this gives 

= A3 - - / dx^\xf^{xf. (IV.IO) 



One can see that -B3 does not commute with Aq, Ai and A2. There were quantum corrections 
even for Ai and ^2, exhibited in the replacement of {H^Y ^i^d {H^Y by Hq{Hq — 1) and 
Hq[Hq — 1){Hq — 2). The corrections are expressed in terms of an operator which was already 
a generator of the sub-algebra of commuting operators. The same treatment of i?3 has given 
an operator which differs from H^^ by a non-commuting part because one of the quantum 
corrections does not arise with the correct coefficient. One might observe that it corresponds 
to a delta function interaction: that is, in Fock space, its effect is only felt at the boundaries 
Xi = Xj. This does not explain the term away, but the explanation is closely related as 
we shall see in the next section. It is useful also to consider the quantum expansion of 
log[exp(— iAL/2)r(A)]. It may be obtained by taking the logarithm of 01.12^ . or by direct 
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computation, which will be given below. Either way the result is 



log(f 



-iAL/2, 



(A)) 



IC 



IC 



3c 



-[H, + cH,~-N\+-\H, + -H, + c^H,--N 



0{\-^). (IV.ll) 



This is not the normal ordered form of the classical expansion f irvTSj) ! Whilst it is true that 
each of the coefficients obtained by the normal ordering recipe is a conserved quantity, up 
to H^, there are quantum corrections to the asymptotic expansion beginning even with the 
Hi term. Moreover, the expansion ( 1IV.11I) is equivalent to (11.121) by simple exponentiation, 
so the normal order of the terms in fllV.lip implies the lack of normal order for ^3. 

The best way to control this ordering problem is to use lattice regularisation. That is, 
we solve the QNLS model on a lattice exactly using the QISM and this will allow us to 



18| for the 



calculate the quantum corrections directly. To do this we need the results of ref. 
inverse scattering scheme for both the classical and quantum cases. The L operator for the 
lattice version has the form 



L{n\\) 



1 - zAA/2 + (cAV2)^t^n ~i^.fci,ip,, 

+il^^cpnilJn 1 + A/2 + (cA V2)^t 



Pr. 



(cAV4)4^„. 



(IV.12) 



Here A is a length (the step length) which we use to take the continuous limit. Also we are 
using a canonical Bose field ipn on the lattice. For the classical model we have the Poisson 
brackets 

= r^~H.mn, {^m^^n} = {^L^l} = 0, (IV.13) 

where 5mn is the Kronecker delta symbol. For the quantum case the equivalent commutation 
relations are 

[V;^, V^t] = A-i5„,„, [^^, = [,1^1, V^t] = 0, (IV.14) 
The transition matrix T(m, ?7,|A) is given by the usual formula for the QISM: 



T{m,n\\) = L{m\X)L{m - 1|A) ■ ■ ■ L{n + l\l)L{n\X), (m > n). 



(IV.15) 



The monodromy matrix r(A) is well defined for a one-dimensional lattice of sites: it is 



T(A)=T(A^,1|A). 
15 



(IV.16) 



The transfer matrix which incorporates periodic boundary conditions is simply the trace: 



r(A) =trT(A). 



(IV. 17) 



The central identity is that 



{r(A),r(/i)} = 0, (classical), 
[r(A),r(/i)] = 0, (quantum). 



(IV.18) 



It is important that it is still valid in the continuous limit. In fact the classical r- matrix, 
which is 

■"l o" 



r(A,/i) 



and the quantum i?-matrix, which is 



\ — 



10 
10 
1 



(IV.19) 



/(/i,A) 








/(/i,A) 

gifx,\) 1 
1 gifx,\) 

/(/i,A) 

/i — A + 



^(/^, A) 



tc 



(IV.20) 



(IV.21) 



fi — X fi — X 

do not depend on A and this is a simplifying factor in taking the limit A — )■ 0. We see 
immediately from these formulae that normal ordering is lost on the lattice. That is. 



T(A) ^ :T^'(A):, r(A) ^ :r='(A): 



(IV.22) 



The crucial point is that exact integrability remains in the continuous limit A — 0, and 
expansions around A = are easy to get. Moreover, in any such expansions the first few 
terms will be normally ordered. This is exactly the behaviour noted above. 

Now we proceed to the calculation for the continuous case in a finite box. The eigenvalues 
6{X) of the transfer matrix r(A) = A{X) + -D(A) are known from the QISM, they are 



N 

~i\L/2 J-j- / ^ 



IC 



X — 



N 



IC 



X — k. 



(IV.23) 
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Here the momenta kj must satisfy the Bethe Ansatz system of equations: 



AkiL 



N 

n 



ki — kj + ic 
ki %c 

1 ' J 



(IV. 24) 



Decomposition of 9{\) in the A — )■ —zoo hmit may now be made by expanding the (finite) 
product of eigenvalues (]IV.23I) . This is similar to eq. (9.3.27) of ref. 15|, here we correct 
two numerical errors in the coefficients. The result is 



-i\L/2 



9{X) 



- A— ioo A A 

N 



N 



- ^ I E ^? + - 1) E - TT^^^ - - 2) 



ic 
A4 



AT 



TV 



TV 



(AT _ i)(iv - 2) 5^ A;, + ^iV(iV - l)(iV - 2)(iV - 3) ) . (IV.25) 



In section Hn] we already constructed operators with these eigenvalues, on a complete set of 
states. This identification leads to the following decomposition of r(A) itself: 



-i\L/2 



r(A) 



1 + A-^o + A-'Ai + \-^A2 + A-^Ag + 0(A-^), 



A— ioo 

where carries the required quantum correction, viz 



(IV.26) 



^3 =:[<^Hf--{H'i 



c\\2 _ ^2uc\uc\ 



+ ^ J dx^^xf^ixf. (IV.27) 



This is the result quoted in eq. f ll.lSp . It follows from the commutativity of the family r(A). 
A similar decomposition may be made for the logarithm of the eigenvalues. Corresponding 
to (1IV.25P we have 



log(e-*^^/2^(A)) 



A— ioo 



TV . \ . / N TV 2 

i-x^'-5(5:*.+>)-S(e*i+-e%-t^ 



TV 



N 



2 1 A3 

N 



,3 



^ ( E *l + ^ E ^1 - E - ) + ■ ('^.28) 
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from which we obtain (lIV.lip . 

We conclude by mentioning that the calculations can also be performed using the methods 



of ref. 



18| . in which a staggered lattice model is introduced in order to make A = —2?/ A a 



special point where the local transition operators become one-dimensional projectors. The 
higher Hamiltonians may then be extracted directly using logarithmic differentiation of the 
transfer matrix at this point. On the lattice, the failure of normal ordering is no surprise 
since as we already noted in eq. (]IV.22p . The calculations are very long and will not be 
given here. 



V. DIFFICULTIES WITH DIRECT ASYMPTOTIC EXPANSION 

We have demonstrated in the previous three sections that when due correction of errors 
is made to formulae given in refs. Q, [l5|, there is no problem with the generation of 
higher conservation laws for the QNLS using either the differential equation formulation or 
the QISM. We now address the question of what goes wrong with the direct asymptotic 
expansion of the operator A{X) for the continuous QNLS in an infinite box. First we give 
some definitions and make some general observations. For a classical field theory involving 
the field ipi^x, t), and for given functions am„(xi, . . . , Xm', Hi, ■ ■ ■ , Vn), we define a functional 
A'^^{ip,ip) of the form 

A''\'ip,ip) = J d'^xd''yarn.n{xi,.-.,x^;yi,...,yn)^p{xi)---i){xrn)'4'{yi)---'^{yn)- (V.l) 

mn 

In a quantum theory, involving the field \E'(x, t), operators A may be constructed similarly. 
In second quantised form we write 

A = J2 d'^xd''ya„,n{xi,...,x^]yi,...,yn)'^\xi)---^\x^)<i/{yi)---<i'{yn). (V.2) 

mn 

Because we must specify the ordering of operators, there are many possibilities for A. Here 
we have shown the "normal ordered" form: we write A = ■.A'^^: to indicate normal ordering. 

There seems to be a folk theorem which says that, whenever we have a Poisson bracket 
relation for classical observables (for instance, a conservation law) then the corresponding 
quantum version must use the normal ordered form of the classical functional. Such a 
connection is not a necessary ingredient for exact integrability of a quantum theory. It is 
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well known 



13| that the classical and quantum coefficients A'^^{X), B'^^{X), A{X) and B{X) 



of the Zakharov-Shabat scheme for the NLS, in an infinite box, satisfy 

AiX) = ■.A^'iX):, (Im(A)<0), 

B\X) = ■.B'\X):, (Im(A) = 0). (V.3) 

For the CNLS in an infinite box it is the expansion of the logarithm of ^'^'(A) which generates 
the higher Hamiltonians in a simple (linear) way. Since the logarithm is non-linear, we would 
expect quantum corrections in this expansion. Thus it causes no difficulty for the QISM 
that we should have 

H2 = -.Hf:, H, ^ -.Hf: . (V.4) 

What is important is that Hf can be recovered from in the quasiclassical limit, and this is 
so because we have shown that there is no discrepancy in the various terms of a^n- Again, 
the expansion of A{X) itself leads to products of quantities in the higher coefficients 
Af^: the difference between normal ordering the if^^ or the Af^ once more involves quantum 
corrections. 

In ref. 15| some problems are indicated with the expansion of A{X). We have already 



mentioned that some of these are computational errors, and we have given the corrected 



formula for A^ in II. 161 The substantial argument given in ref. 15| is that, if the asymptotic 
methods used for the decomposition of /l'^'(A) are repeated with ■.A'^^{X): with the normal 
ordering retained at each step, then we should get the expansion coefficients as An = :A'jl: 
for all n. This is not so, and we need to see why the analysis fails for n > 3. We stress 
that the manipulations used in the quoted analysis are purely formal, and the calculation 
of the quantum corrections (which appear from the normal ordering) depends on using the 
canonical commutation relations for fields \&^(x) and "^{y) in integrals which have x = y as 
one limit. While this kind of formal analysis may work well in many cases, we have no right 
to expect this. Any proper asymptotic analysis will depend on the action of A{X) in Fock 
space as an integral operation. In the QISM, A{X) is formally defined by its formula in the 
second quantised form jl2|: viz 

A{X) = ^c" dx^'dy"' 6{xi < yi < ■ ■ ■ < xn < yN) 

n>0 

exp [iX{xi -yi + --- + XN- 2/7v)] ^K^i) ■ ■ ■ -^\xN)-^{yi) ■ ■ ■ ^(y^v)- (V.5) 
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Here 6{xi < yi < ■ • ■ < xn < Vn) stands for the indicator function of the set |xi < yi < 
... < xn < Vn}- The action of this operator in Fock space was given in ref. 13|. Let 
I/) and \g) be two N particle states specified by symmetric functions f{xi, . . . ,xn) and 
g{xi, ...,xn) via 

|/) = j dx''f{xu...,XN)'^\x,)---^\xNm, (V.6) 
with a similar equation for l^f), then the action 

\9)=Mm) (V.7) 
is given by the following integral operator 

g{xi < ■ ■ < Xn) = fixi, . . . , Xn) + XI XI / ■ ■ ■ d$,i 

n=l h<-<in n 

exp[iA(xi, -^iH hXi„ - Cn)] /(^i, • • • , Cn|6 Xi^,. . . ,in^ Xi^). (V.8) 

Here the notation for the integrand means that the indicated changes of variables are made in 
the function /. Also the evaluation of g{xi, . . . , xn) for orderings other than xi < ■ ■ ■ < xn 
is by symmetrisation. 

An integral operator typically represents a boundary value problem. Direct computation 
from flV.81) shows that this is so for A{X). The functions f{xi, . . . , xn) and g{xi, . . . , xn) 
are related by the boundary value problem 

n( + ^fH l^i^^^ ...,xn) = YI[\ + i— - ic]fixi, . . .,xn), (V.9) 



eg + 



d 



dxj 
d 



^dxj dxj^i 
It is shown in ref. 







C/+ ( 









dXj 

d 



d 



dxn dx 



/ 



(V.IO) 



- a;,-+i=x,-|-0 



15| that the various operators of the QNLS theory (c 7^ 0) are intertwinings 
of corresponding free operators (c = 0) restricted to a domain in which appropriate boundary 
conditions are satisfied. Eq. fIV.lOp tells us that A{\) may be restricted to the appropriate 
domain: that is, A{X] preserves just the correct boundary conditions. In fact we can see 



this directly from ref. 13| where it is shown that the integral operator fIV.Sp is diagonal on 
the Bethe-Ansatz eigenstates. The latter are a complete set among precisely those functions 
which satisfy the boundary conditions required to define the commuting operators Hn for 
the interacting case as restrictions of the free Hamiltonians. 
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The decomposition of the operator (IV. 80 in inverse powers of A, using the usual techniques 
of integration by parts, gives a non-uniform asymptotic expansion in the variables xi, . . . ,xi\i, 
which fails exactly at the boundaries Xj = Xk- The two particle sector will suffice to illustrate 
the ideas and in fact we only need go to 1/ to see how things work: 

9{x,y) ^ f{x,y) 

/ 2 1 1 

+ [{U^. y) + /.(^, y)] - e^'^^-^H y) + 2/)}] ) + ■ ■ • • (v.n) 

Away from the boundaries x = this expansion correctly gives the differential parts of the 
operators in the asymptotic expansion, since we may neglect the exponentially small correc- 
tions when X < y and A — )■ —zoo. The results then are the same as for the non-interacting 
theory. To complete the expansion we must find out what happens at the boundaries and 
this cannot be deduced from (IV. lip . However, our comments above show us that we may 



identify the operators in the expansion of A{\) by their differential parts found from the 
asymptotic expansion away from the boundaries Xi = xj. These give us the free Hamiltoni- 
ans: the intertwining property takes them into the corresponding interacting Hamiltonians. 
The details of this calculation are equivalent to the calculations involving eigenvalues given 
in section HVl so we do not repeat them here. The result is therefore that the correct asymp- 
totic expansion of A{X) in terms of higher Hamiltonians is given by eqs. (II.12l) - (II.16p . in 
agreement with the calculations made from a lattice limit. 



VI. CONCLUSIONS 



As we mentioned in the introduction, there have been a number of papers which have 
raised various mathematical questions about the QISM solution of the QNLS. In this paper 
we have been concerned with the most serious objections, which suggested that the conser- 
vation laws are flawed. We have shown that they are not. However, it must be stressed 
that explicit formulas for the higher conserved quantities are difficult to get and to use be- 
cause one must go through singular calculations. While these difficulties may impair their 
practical utility, it is certainly not a flaw in the QISM, and that is the chief concem of this 
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paper. Fortunately, there exists a well behaved lattice regularisation of the model which can 



control these problems. The same comments apply to the quantum trace identities |2l|. In 
this view, everything in the continuous case is understood as the appropriate limit from the 
lattice. This controls the ordering problem for these laws, and shows that normal ordering 
is not correct beyond H^: also that there exist quantum corrections beginning with As in 
the expansion of A{X). 

Our conclusion is that the Bethe Ansatz solution and the QISM give the same (valid) 
conservation laws: also that the quasi-classical limit is correct. So there is no failure of 
the QISM for the quantum non-linear Schrodinger equation. This was a most important 
point to resolve now that the QISM seems poised to solve the long-standing problem of the 



construction of correlation functions for solvable models 
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